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ABSTRACT 

Results from the first two years of data from the Taiwanese- American Occultation Survey (TAOS) 
are presented. Stars have been monitored photometrically at 4 Hz or 5 Hz to search for occultations by 
small (~3 km) Kuiper Belt Objects (KBOs). No statistically significant events were found, allowing 
us to present an upper bound to the size distribution of KBOs with diameters 0.5 km < D < 28 km. 
Subject headings: Kuiper Belt, occultations, solar system: formation 



1. INTRODUCTION 

The study of the Kuiper Belt has exploded si nce the 
discovery of 1992 QB1 by Uewitt fe Luul (j!993f) . The 
brightness distribution of objects with R magnitude 
brighter than ~26 is relati yely well-esta b lished by many 
surveys, most recently by iFraser et al.l (|2008t and ref- 
erences therein) . The brightness distribution is ad- 
equately described by a simple cumulative luminosity 
function £(< R) = IQO'(R-Ro) deg~ 2 , where R ~ 23 and 
a ~ 0.6, for objects with magnitude R < 26. There is 
clear evidence for a break to a shallower slope for fainter 
objects: the deepest survey, conducted using the Ad- 
vanced Cainera 1 fpj_Surveys on the Hubble Space Tele- 
scope (|Bernstein et al.ll2004h extended to R = 28.5, and 
found a factor of ~25 fewer objects than would be ex- 
pected if the same distribution extended into this range. 

The size distribution of Kuiper Belt Objects (KBOs) 
is believed to reflect a history of agglomeration during 
the planetary formation epoch, when relative velocities 
between particles were low and collisions typically re- 
sulted in particles sticking together, followed by destruc- 
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tive collisions when the relative velocities were increased 
by dynamica l processes after the gian t planets formed 
(ISternl fl99l iDavis fc Farinellal fl997t IStern fe Colwelll 
19971: IKenvon fc Luulll999aU bl: iKenvon fc Bromlevl 12004 
Pan fc Saril l2005h " The slope of the distribution func- 



tion for larger objects reflects the early phase of ag- 
glomeration, while the shallower distribution for smaller 
objects reflects a subsequent phase of destructive col- 
lisions. The location of the break moves to larger 
sizes with time, while the distribution for smaller ob- 
jects is expected to evolve towards a steady state col- 
lisional cascade (|Kenvon fc Bromlevl 12004 iPan fc Sari 
I2005T ) . Models for the sp ectrum of small bodies differ 
between fPan fc Saril (l2005f). who derived a doubl e power- 
law distribution, and IKenvon fc Bromlevl (|2004) . whose 
simulations show more structure, depending on material 
properties. 

Thus, the size spectrum encodes information about the 
history of planet formation and dynamics. However, the 
size spectrum for small KBOs is not constrained by the 
imaging surveys because the objects of interest are too 
faint for direct detection using presently available instru- 
ments. These small objects may, however, be detected 
indirectly when they pass between an observer and a dis 
tant star (lBailevlll976tlDvsonlll99llAxehod et all 
Brown fc Websterl Il997t iRoaues fc Moncuauetl 



1992 



2000: 



Coorav fc Farmeril2003l : lNihei et al.ll2007l ). The challenge 



confronting any survey exploiting this technique is the 
combination of very low anticipated event rate and short 
duration of the events (typically < 1 second). 

Other groups are a t tempt ing similar occultation sur- 
veys. IRoaues et al.l (|2006l) reported three events in 
10 star-hours of photometric data sampled at 45 Hz, 
which they modeled as objects at 15 AU, 140 AU, and 
210 AU, respecti vely, placing the i nferre d objects outside 
the Kuiper Belt. iBickerton et alj (|2008l ) reported results 
of 5 star-hours of data sa mpled at 40 Hz, dur ing which 
no events were detected. IChang et al.l {2006) reported 
a surprisingly high rate of possible occultation events in 
RXTE x-ray observations of Sco X-l, but many of these 
events have since been attributed to i nstrumental effects 
(|Jones et alj|2008t IChang et al.ll2007ft . 

We report here the first results of the Taiwanese Amer- 
ican Occultation Survey (TAOS). TAOS differs from the 
previously reported projects primarily in the extent of 
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Fig. 1. — Demonstration of detrcnding by the filtering process. The left panel shows a raw lightcurve set f(t), and the right panel shows 
the lightcurve set h(t) after filtering. The total number points in the lightcurve set is 26,622, corresponding to 89 minutes. Here only one 
of every 10 points is plotted for readability. 

the photometric time series, a total of 1.53 x 10 5 star- 
hours, and in that data are collected simultaneously with 
three telescopes. Some compromises have been made in 
regard to signal-to- noise (SNR), which is typically lower 
than in previously reported surveys, and in cadence, 
which is 4 Hz or 5 Hz, in contrast to the higher rates 
mentioned above. The substantial increase in exposure 
more than compensates for the lower cadence and SNR, 
and we are able to probe significant ranges of the model 
space for small KBOs. We have also developed a statis- 
tical analysis technique which allows efficient use of the 
multi-telescope data to detect brief occultation events 
that would be statistically insignificant if observed with 
only one telescope. 
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2. DATA AND ANALYSIS 

TAOS has been collecting scientific data since 2005. 
Observations are normally carried out simultaneously 
with three 50 cm telescopes (A, B, and D, separated 
by distances of 6 meters and 60 meters; this system is 
described bv lLehner et~aTl (|2008af l). Over 15 TB of raw 
images have been taken. We report here on the first two 
years of data taken simultaneously with all three tele- 
scopes. The data set comprises 156 data runs, where a 
data run is defined as a series of three-telescope observa- 
tions of a given field for durations of ^90 minutes. Thirty 
data runs with a 4 Hz sampling rate were taken before 
2005 December 15, and 126 data runs were collected sub- 
sequently with a sampling rate of 5 Hz. Only fields with 
ecliptic latitudes \b\ < 10° were analyzed. Over 93% 
of the data were collected in fields with |6| < 3°, so 
the results of our analysis are relevant to the sum of 
the c old and excited KBO populations (|Bernstein et al.l 
120041 ). No data run was included unless each star was 
sampled more than 10,000 times in each telescope. The 
angle from opposition in these data runs is distributed 
from 0° to 90°. The number of stars (with R < 13.5, 
which typically gives a SNR > 5) monitored in the data 
runs ranges between 200 and 2000 1 . 

The images were analyzed usi ng an aperture photom- 
etry package (Z hang et alj 12008ft devised exclusively for 
TAOS images. Lightcurves were produced for each star 
by assembling the photometric information into a time 
series. A star in each data run has a lightcurve from each 
of the three telescopes. The data presented in this paper 
comprises 110,895 lightcurve sets (where a lightcurve set 

Information on the TAOS fields i s available at 
http : //taos . asiaa.sinica.edu.tw/taosfield/ 



Fig. 2. — Rank-rank plots of the lightcurve set shown in Fig.[T] 
Each point shows the ranks at a single time point for telescopes A 
and B. The ranks from the raw lightcurves (Fig.[T] left) are shown 
in the left panel, and the ranks from the filtered lightcurves (Fig.[T] 
right) are shown in the right panel. 

is defined as a set of three lightcurves, one for each tele- 
scope, for the same star in a single data run) , containing 
7.1 x 10 9 individual photometric measurements. 

The photometric data are not calibrated to a stan- 
dard system. Changes in atmospheric transparency dur- 
ing a data run produce flux variations (Fig.[TJ) that could 
undermine our occultation search algorithm. Such low- 
frequency trends in a lightcurve can be removed by a 
numerical high-pass filter that preserves information of 
a brief occultation event, typically with a duration less 
than 1-2 data points with the TAOS sampling rate. Our 
filter takes a time series of fa measured at_time ti to 
produce an intermediate series gi = fa — fa where fa 
is the running average of 33 data points centered on 
ti. The series Qi is then scaled by the local fluctuation, 
hi = 9il<y{gi) with <r(<fc) being the standard deviation 
of gi of 151 data points centered at U. Both the mean 
and standard deviation are calculated using three-sigma 
clipping. This filtering proves effective to remove slow- 
varying trends in the lightcurve, while preserving high- 
frequency fluctuations that we aim to detect, as illus- 
trated in Fig. [I] 

We now confront the two central challenges in the 
search for extremely rare occultation events in these 
data: (1) to search for events simultaneously in three 
parallel data streams, and (2) to determine the statisti- 
cal significance of any rare events that are found. The 
second is not straightforward because the statistical dis- 
tribution of our photometric measurements is not known 
in advance; approximations based on Gaussian statistics 
are unreliable far from the mean. This motivates a non- 
parametric approach. 

We thus found it useful to represent each data point 
by its rank in the filtered, rescaled lightcurve data hi. 
That is, the rank of a data point ranges from r = 1 (low- 
est h) to r = Ap (highest h), for a data run comprising 
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Fig. 3. — Histogram of r\ for a lightcurve set implanted with a 
simulated occultation by a 3 km KBO. Each lightcurve comprises 
26,637 points and the ranks from the three telescopes are found to 
be 1, 3, and 1. The event is clearly visible at r\ = 29.47, and the 
probability of a rank product of 3 or lower from random chance 
is P = 3.7 X 10~ 13 . The overlapped dash line is the theoretical 
distribution of r\. The theoretical and implanted lightcurves are 
shown in the insert, offset vertically for clarity. 

Np photometric measurements taken with telescope A, 
B, or D. The rank triples (rf ,rf ,rf) form the basis of 
further analysis of the multi-telescope data. The statis- 
tical distribution of these ranks is known exactly, since 
each rank must occur exactly once in each time series. 
Thus, the probability that a given rank will occur at 
time U is P(ri) = l/Np. When the photometric data 
are uncorrelated, the probability that a particular rank 



triple will occur is simply P(r^ 



„B n 



) = l/Np*. This 



allows a straightforward test for correlation between the 
photometric data taken in the three telescopes: the rank 
triples should be distributed uniformly in a cube with 
sides of length Np. A non- uniform pattern in the cube, 
on the other hand, indicates correlation. Fig. [2] shows 
an example of the rank series of one telescope against 
another; this indicates that the raw photometric data fi 
are strongly correlated, but the filtered data hi are not. 

Given that the rank triples are uncorrelated, the ranks 
can be used to search for possible occultation events, as 
follows: A true occultation event will exhibit anomalous, 
correlated low ranks in all three telescopes. The rank 
triples thus allow an elegant test for the statistical sig- 
nificance of a candidate event. Consider the quantity 
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Since the ranks are uncorrelated (unless we have an oc- 
cultation event), we can calculate the exact probability 
density function 2 for r\. This in turn allows us to com- 
pute the probability that a given triple of low ranks ran- 
domly occurred in an uncorrelated lightcurve set; this is 
our measure of statistical significance. An illustration of 
the power of this approach is shown in Fig.[3j where a 
simulated occultation event is readily recovered. 3 

We thus screened all of our series of rank triples for 
events with low ranks in all three telescopes. In the anal- 
ysis reported here, we considered only events for which 
P(v > c) < 10~ 10 , which leads to an expected value 
of 0.24 false positive events in the entire data set of 
2.4 x 10 9 rank triples. No statistically significant events 

For small n, this distribution can be approximated by a V dis- 
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-''/(JVt-I)!, where 2V T 



tribution of the form P(rj) 
is the number of telescopes. 

3 Details of our statistical methodology will be described in 
ILehner et all Il2008bl) . 
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Fig. 4. — Effective solid angle of the TAOS survey as a function of 
KBO diameter D (solid line), and the product of the upper limit of 
the differential surface density by the effective solid angle (dashed 
line, in arbitrary units). 

emerged from this analysis. 4 

3. EFFICIENCY TEST AND EVENT RATE 

While the example event shown in Fig. [3] is readily 
recovered, objects with smaller diameters or which do 
not directly cross the line of sight might not be so eas- 
ily detected. An efficiency calculation is thus neces- 
sary for understanding the detection sensitivity of our 
data and analysis pipeline to different event parame- 
ters, notably the KBO size distribution. Our efficiency 
test started with implanting synthetic events into ob- 
served lightcurves, with the original noise. Each origi- 
nal lightcurve fi was modified by the implanted occul- 
tation, ki = fi — (1 — di)fi, where di is the simulated 
event lightc urve (w ith baseline di — > 1 far from the event, 
(jNihei et alj|2007l )). and is the average of the original 
series over a 33 point rolling window. Since we preserved 
the original noise in the modified lightcurves, the noise 
where the implanted occultation event takes place — for 
which the flux diminishes — would be slightly overesti- 
mated, hence our efficiency estimate is conservative. 

We assumed spherical KBOs at a fixed geocentric dis- 
tance of A = 43 AU. (Given our sampling rate, vary- 
ing the KBO distance within the Kuiper Belt has little 
effect on our simulated lightcurves.) The event epoch 
to was chosen randomly and uniformly within the dura- 
tion E of the lightcurve set. The angular size 6** of each 
star, necessary for the simulated lightcurve calculation, 
was estimated using stellar col or and apparent m agni- 
tude tak en from the USNO-B (|Monet et al.ll2003h and 
2MASS (jSkrutskie et al.ll2006h catalogs. The impact pa- 
rameter of each event was chosen, again randomly and 
uniform ly, between and H /2, where H is the event cross 
section (jNihei et al.ll2007j ). 



H = 2 [(\/3>)f + (D/2)i 



.A. 



Here D is the diameter of the occulting object, and F 
is the Fresnel scale, F — y AA/2, where A = 600 run is 
the median wavelength in the TAOS filter. The relative 
velocity between the Earth and KBO u re i, necessary for 
the conversion of the occultation diffraction pattern to 
a temporally sampled lightcurve, is calculated based on 
the angle from opposition during each data run. 



4 A candidate event was reported in Chen et al. (2007). This 
event had a significance of 3.7 X 10 — 10 , which did not pass our cut 
on rj. We expect to have ~1 false positive event at that significance 
level or higher. 
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Fig. 5. — TAOS upper limit to the size spectrum (solid line) 
assuming a power-law distribution; double power-law fit from 



IBernstein et alj (2004) (dashed line, e xtrapolation at D < 28 km 
shown in dotted line). Results from Bernstein ct al. (2004) are 
shown as d ata points. Upper limits from IBickerton" et al.l (120081 1 
(BKW) and I Jones et~aH (2008T) (RXTE) also shown. An albedo of 
4% is assumed for computing magnitude. 

To adequately cover a wide range of parameter space, 
two efficiency runs were completed in which we implanted 
each lightcurve set with exactly one simulated occulta- 
tion event. For each event, the diameter of the KBO was 
chosen randomly according to a probability, or weight- 
ing factor, wd- In the first run, objects of diameters 
D = 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9 km were added with 
wd = 1/6 for each diameter. In the second run, ob- 
jects of diameters D = 1, 2, 3, 5, 10 and 30 km were 
added with weights w D = {100, 100, 100, 30, 5, l}/336. 
The modified lightcurves ki were reprocessed using the 
same procedure described in The recovered events 
and the event parameters were then used to calculate 
the number of expected occultation event in our survey. 
That is, we calculated the quantity 



n a (D) 



the model- dependent size distribution dn/dD, and the 
model-independent effective solid angle Q C (D), which de- 
scribes the sensitivity of the survey to objects of diam- 
eter D. Given the observed number of events and the 
value of il c (D) resulting from the efficiency calculation, 
we can place model-dependent limits on the the popula- 
tion of KBOs. Based on the absence of detections in this 
data set, any model with a size distribution such that 
N cxp > 3.0 is inconsistent with our data at the 95% 
confidence level. 

Note that there are an infinite number of models that 
satisfy the above requirement. We thus make the rea- 
sonable choice of a power-law size distribution dn/dD — 
Ub(D/28 km) -9 , where «b is chosen such that the cumu- 
lative si ze distribution is contin uous at 28 km with the re- 
sults of IBernstein et ail (|2004D . We integrate Eq.[T]from 
£>2 = 28 km down to our detection limit of D\ — 0.5 km, 
and solve Eq.[T]with iV GX p = 3, to find q = 4.60. Our null 
detection thus eliminates any power law size distribution 
with q > 4.60 at the 95% c.L, setting a stringent upper 
limit (see Fig.[S]) to the number density of KBOs. 



4. CONCLUSION 

We have surveyed the sky for occultations by small 
KBOs using the three telescope TAOS system. We have 
demonstrated that a dedicated occultation survey us- 
ing an array of small telescopes, an innovative statistical 
analysis of multi-telescope data, and a large number of 
star-hours, can be used as a powerful probe of small ob- 
jects in the Kuiper Belt, and we are thus able to place the 
strongest upper bound to date on the number of KBOs 
with 0.5 km < D < 28 km. We continue to operate 
TAOS, soon with an additional telescope, and will re- 
port more sensitive survey results in the future. 



where the sum is over all lightcurve sets where a sim- 
ulated event is successfully recovered in the reanalysis 
(Fig. 2]). Essentially £l c (D) is the effective solid angle of 
our survey, insofar as TAOS can be considered equiva- 
lent to a survey that is capable of counting every KBO 
of diameter D in a solid angle fl e with 100% efficiency. 

The expected number of detected events by KBOs with 
sizes ranging from D\ to Di can then be written as 



(i) 



where dn/dD is the differential surface number density 
of KBOs. The integrand of Eq.Q] contains two factors: 
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